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ABSTRACT 


The military threat of hostile submarines is increasing and the need for effective 
Anti-Submarine Warfare (ASW) operations is also increasing. In response, the ROK 
government and military have improved their ASW capabilities. In this thesis, the 
recommended track spacing for an Archimedes spiral search in a datum search problem 
was studied. To find a recommended track spacing, three analytical approaches were 
explored. Each of three analytical approaches has its own strengths and weaknesses. This 
analysis combined three analytical functions into a single parameterized expression. To 
find the best-fit parameters maximizing probability of detection, a simulation experiment 


with a NOLH (Nearly Orthogonal Latin Hypercube) design was used. 


THIS PAGE INTENTIONALLY LEFT BLANK 


vi 


TABLE OF CONTENTS 


I. TNTRODUC TION vivcceccviscecsteckevactcasieaucdecavedeudacectaesduenevbalcssevacceseuvedeousessiéesavecsverexdicoies 1 
A. BACKGROUND irieseciesscavccarievssidiectuccssvesstacseiesissytnersccsmmivacvaaacsaensipiawadae 1 
B. SCTNANCIO) scsi visicesaccdisnshceieiarasaiets eedticteses aeuateniied aedaranebds teed eae avielaek 1 
C. RESEARCH QUESTION scsiciscicsdiscesistsiccvscdeiecsetisivesstcosconstécasvsecdaveebocddxnesteees 2 
D. STRUCTURE, OB“ THESIS ssisdicsndévesits vacceaceeesscscceccedsvievesdeciescotstts vecessdesceiseiesia 2 
II. SEARCHER MOTION MODEL vicccccccsccossesscousscssseneesscevvesoncsdetesceonsncsatesenssevenssesnseses 3 
A ARCHIMEDES SPIRAL SEA RCH1,..........ccccccssssscsssccsscccssscccesscccssscceesscsesees 3 
TH: TARGET MOTION MODEL svscscscsseccscsscecsascacsenscscincscsconsestestsesdeieseasvsstsdasesscaesevetecas 5 
A. RANDOM TOUR besesceadisscdasccccasvadencdscacsceusstevsidsdccebevsatsvsedesieecssialerindsteceusvtcne 5 
B. PROBEEM PA RAME TERS sivicscscavavsessavvnsdesssowvucvedvsvaseereavdes cavevinsituevaccudisevucs 5 
1. UCR ar Bet Speed ) ey cssecevaces cxcicdvv eres cess cacuatetetacousenicaneevaveciv dae hieesbcasoaves 5 
2. A. (Course Change Rate). sscsisscivesecdvaccpecesasceaseasvaredaiveansnecenocteonaveccnesye 6 
3. GC (ODSETVAtON ELLOL) seccisscsstscsecesessosescecesecssossscossecsecccsossosesescscousbesssee 6 
C. NORMAL APPROXIMATION FOR THE PURE RANDOM TOUR 
TARGET ieaieteticincaiat in auainntacinienniias einen 6 
IV. ANALYTICAL APPROACHS TO FIND GOOD TRACK SPACINGS ............. 11 
A. IT aaa si TRY faa Pe A RY RE PON ORES ROAR NAGY IPRA EPR RPI MOR 11 
B. Sipay taesseaseesesccssecsecessscnscessscnssssssensssssssessecessnenecsssscesessussessssassessesessnsnecessenecosencecs 13 
C. SDI aap teeissiasavaabucsiviashcaivsapheueleatinheealehetestatoesetaniee ice biecsaGaaeeianas 16 
D. PROPOSED EXPRESSION ssicscscsscssnssdisserscuscvescisdessctcadcassencaasestcdacnscdcinessees 16 
V. DESIGN OF EX PEREVIEN Fiivevessvesseacsvedas saseasecusscdusatetesesedtvasdescvacded carves decenselavistvrs 17 
A. NOLH (NEARLY ORTHOGONAL LATIN HYPERCUBE) DESIGN...17 
B. THE MAXIMUM PROBABILITY OF DETECTION ( Pd,,,,) FOR 
BACH: DESIGN. POIN T sstsssisssesistsces caus sasagasiasencensconcauedosusccvsceneesseiovecuniteseavune 18 
C. FINDING THE BEST G8, AND sisccviesi sieciavicdevsdssetentdbidiicseseivecdeiedecetes 19 
Wie. “SIMULATION RESULTS sis seisciecncaeteisiarsrdncasivdcraratealicontaeedacienianic 21 
A. IWEAUIN: RES UIE Dvsssicspsnutitcsvicavacevcicesdasacaicuecsuasaduastesassescsadseeetacassascesi vecutatersies 21 
B. VISUALIZATION OF MAIN RESULT USING CONTOUR PLOTS.....21 
VII. CONCLUSION / FUTURE RESEARCBH.uuu.........ccceccssscecescccesscccesccccesccecssccecsscceceese 25 
A. CONE USION ss issasavecaivcaubdesnsnskcunteachcdsvse ceases onndiats cesdeanicueh cdeatviecsiviusbicierstesate 25 
B. FUPU RE RES EA RC Ei eoiscesetsscatsicie vectors cates shasta .cenatets catanasaneteaeiensiate 25 
APPENDIX A. NOLH (NEARLY ORTHOGONAL LATIN HYPERCUBE) 
DESIGN FOR FIVE BLOCKING FACTORS ...0.........ccssssscssssscssscccesscccsssccesscccenees zt 
APPENDIX B. JAVA CODE FOR DATA COLLECT ION .............ccsssccssssscssreccssccceessceees 29 
LIST OF REFERENCES wisctaiacsssccisasscasdssackiandececeaestacadinedcleaecssecksivescteoseevevataicseseaneepecdaiacocies 31 
INITIAL, DISTRIBUTION LIST ssssecssssccssccesessedebscdsccescensiteerscsssevsosssiessisesseataiaisssessdsceesedsoiss 33 


THIS PAGE INTENTIONALLY LEFT BLANK 


Vill 


Figure 1. 
Figure 2. 
Figure 3. 
Figure 4. 
Figure 5. 
Figure 6. 


Figure 7. 


Figure 8. 
Figure 9. 


Figure 10. 
Figure 11. 
Figure 12. 
Figure 13. 
Figure 14. 


LIST OF FIGURES 


AS W Sear elt Ste nario. socciyst oc c2ssetacaitaies cconlecaetasten <2 ues secdids ies ateenaeliat acta cancel 1 
RoJers1 £6 AEGIS V0 4cd 11 0 cyanea eer tr ane OPM OT Nant Rr CoM eT en ear tren nano tree 3 
Description of Random Tour Target... eee ceeeeeseecsseceseeeeeeesseecsaeenseeseeeensees 5 
The Initial Rando Pour e 557s see sco' sessed cas ec seaevinensecenetteoes ewes eae 6 
CAlcUlagon Ol 5, sai borsuaaihnnta cemonanah acd aeieteva waned eudloaieed 12) 
Density Function of Radial Position of Target at Time t, Assuming a 

Normal Approximation to the Random Toul.............eeeeeesecsseceseeeseeeeseeenaeenes 14 
KCAlCULAMOMGOE 3,8. ita, wh cunttayanassnscavntie ona tsadatiaated mea sannivgere amauta valet insuas 15 
Crossed DGSiomTOr Sie HACEOIS 5 ca4 5c jh csancsv bead ces scateacuy accsee sacees ooybaseecduateesasbcates 17 
Plot of Estimated Simulation Probabilities of Detection for Blocking 

Factor Design Point. ajcsiissiecsavisceasiedserceasnniaa cs ssedancasashaccons bedava nageadevdesdaativeseadea 18 
Contour Plet.for a “and: 2 Siven: 7 =O. weccssshss ioieciesecreds Wisantoeeetsacs-sabetertees 21 
Contour Plot tor a and: PF S1ven 7H yas sases ves s eosecetequ as easeeevpensuee veer heats, 22 
Contour PIGttor @ ANd: 2 SIVEN HS 2 os ccacdacezetiet ange usasaadenaidaccanentieacheagesandes 22 
Contour Plot for @ and: # Piven 7= 3 i.sasascssergesetaspevcdashieevmensnentns 23 
Contour Plot.for @ and: {given 7 = 4 uiciiccteeiaiici esas 23 


1X 


THIS PAGE INTENTIONALLY LEFT BLANK 


Table 1. 
Table 2. 


LIST OF TABLES 


Probability of Detection ( Pd ) for Each Blocking Factor Design Point 
Table for S,,,, and Pd 


best best 


Xl 


for Each Blocking Factor Design Point......... 


THIS PAGE INTENTIONALLY LEFT BLANK 


xii 


ACKNOWLEDGMENTS 


I would like to express my gratitude to those who contributed to the successful 


and timely completion of this thesis. 


First, I would like to thank my advisor, Professor James Eagle for his patience, 
guidance and constant support. I would like also to thank my second reader, Professor 
Paul Sanchez for his expertise, cooperation and willingness reading this thesis. Other 


members of the NPS faculty also have my appreciation for their excellent instruction. 


My deepest gratitude also goes to the Korea National Defense University, 


especially Professor Jae-young Lee who encouraged and supported me for last one year. 


Finally, I would like to thank my family (parents, sisters, and brother), fiancée, 


and friends for supporting me. 


Xlil 


THIS PAGE INTENTIONALLY LEFT BLANK 


XIV 


I. INTRODUCTION 


A. BACKGROUND 


The Republic of Korea (ROK) is surrounded by several powerful countries and 
also by water on three sides. For several decades, the ROK has faced an increasing threat 
of hostile submarines. In response, the government and military have improved their anti- 
submarine warfare (ASW) capabilities. Weapon systems such as the P-3C maritime 


patrol aircraft (MPA) have been deployed. 


These ASW systems are important for preventing terrorist invasions by submarine 
and to protect the sea lanes. This thesis will develop an ASW search tactic for an MPA 


aircraft searching for a patrolling submarine. 
B. SCENARIO 


The periscope, communications antenna, or snorkel of a target submarine is 
detected by the long-range, surface search ASW radar on an MPA aircraft. The initial 
radar detection provides an estimated target position, called datum. The aircraft then flies 
directly to datum and conducts a spiral-out ASW search for the target. The target 
submarine submerges below periscope depth immediately after radar detection and 


continues its underwater patrol. 


TV. 


transit Vas 
-_- 


ig 


search 





Figure 1. | ASW Search Scenario 


C. RESEARCH QUESTION 


Referring to Figure 1, the problem parameters are: 


U 


V 


(hr ) 


A 


: target speed (kt ) 
: aircraft speed (kt ) 
: range at which the aircraft can detect the submarine (nm) 


: time required for aircraft to fly from where detection occurred to datum 


: the Poisson rate at which the target changes course (hr7') 


The main research question investigated here is as follows: Given problem 


parameters U, V, R, t , and 2, what is the aircraft search track spacing S° which 


maximizes the probability of detection of the submarine by the ASW aircrafts? 


D. STRUCTURE OF THESIS 


Chapters II and III discuss the mathematical model of a target submarine and a 


searcher aircraft which conducts an ASW search operation. Chapter IV develops three 


analytical approaches to suggest good track spacings for the spiral-out search. Chapter V 


develops the design of the simulation experiment and the analysis for constructing a 


meta-model to find the recommended track spacing. Chapter VI discusses the simulation 


results. Chapter VII provides conclusions and recommendations for further research. 


UH. SEARCHER MOTION MODEL 


A. ARCHIMEDES SPIRAL SEARCH 


Aircraft arrives at datum and 


Initial target detection SY Search bepins 
ee t=0 


Figure 2. Search Timeline 


The target is initially detected at time f=—7. At time 0, the searching aircraft 
arrives at datum and begins a spiral-out search with constant track spacing S and speed V. 


The aircraft carries a detector with certain detection range R. 
The aircraft track will follow an Archimedes Spiral [1]. For such a search path, 
the radial distance from datum in terms of the angular position @ in radians is 


OSs 
0) =—. 
ss 20 


So when @ increases by 27 (one revolution), the radial distance increases by S. 


In Chapter IV, it is shown that at time t 20, the approximate radial position at time ¢ of 


a searcher with speed V is 


r(t)= lay 
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Hl, TARGET MOTION MODEL 


A. RANDOM TOUR 


This model assumes the target’s motion is a random tour [2]. The target’s motion 
is assumed to be independent of the search, and not specifically directed toward escape. 
The target may move a considerable distance from the point of initial detection (datum) 
as the search for it proceeds. The target motion can be beneficial for the searcher in that it 
increases relative speed, but it can also be harmful as it increases the size of the area in 


which the target can be located. 





Target Motion Model(U=20, Lambda=0.5) 


29 
29 


60 









































Figure 3. | Description of Random Tour Target 


B. PROBLEM PARAMETERS 


Figure 3 shows an example of a random tour target motion with 4 of .5 course 


changes/hr, and U of 20 nm/hr. 
1. U (Target Speed) 


Target speed U is assumed to be constant. 


2. A (Course Change Rate) 


The target changes its direction at times determined by a Poisson Process [3] with 
rate 7. For instance, if 2 =2 / hour, the target changes its course two times per hour on 
the average, and the time between course changes is exponentially distributed with mean 


1/2 hour. 
3. o (Observation Error) 


Immediately after the initial radar detection at time t=—7, the target’s location 
distribution is circular bivariate normal with equal means (0, 0) and variances (o’,0°) 


In this thesis, we will assume o =1nm, consistent with a long-range radar detection of a 


small target. 


cC. NORMAL APPROXIMATION FOR THE PURE RANDOM TOUR 
TARGET 


For large enough time fr, the Central Limit Theorem [4] guarantees that the target 


position will be approximately circular bivariate normal. 


U : target speed 


T : duration of 1“ leg 
(exponentially distributed with 


mean yy ) 


@ : Target’s course during 1“ leg 
(uniformly distributed between 
0 - 27 radians) 


D : length of 1“ leg 
(exponentially distributed with 


mean U7 ) 





Figure 4. The Initial Random Tour Leg 


To see this, consider the following development. Since target’s course change rate is 1, 
the time to the first course change, 7, is an exponentially distributed random variable 


with rate 1. We write this as T ~ exp(/). So, 


Thus, D ~ exp(4/,). That is D, the distance travelled by the target during the 
first leg, is an exponentially distributed random variable with mean Uy, : 


Now, we want to find the expected value and variance of the target position after 


one leg. Each new course @ is uniformly distributed. So, 
0 ~ U (0, 2r), and 


X =Dcosé. 


Now by symmetry of the cosine function, 


So, 


=£(°) 


Now conditioning on the target course 0, 


E(X*)= f £(x?|0=c) f,(x)ae 


t=0 


And by the symmetry of cos’ (z), 


= (=) feos cE(D*)dr. 


And from the properties of the exponential distribution, 











Now by x, y symmetry, 


Let X,,, be the x coordinate of the target position at time f. 
Then, 
Xoo =X, +X, ++-+Xy, where N ~ Poisson( Ar). 


Since variances add for independent random variables, 





Now, conditioning on N , 


Thus, 


And by the Central Limit Theorem, after a sufficiently large number of turns, we can 


approximate the target’s position with a circular, bivariate normal random variable with 


mean (0,0) and variances in the x and y directions of U Uh 
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IV. ANALYTICAL APPROACHS TO FIND GOOD TRACK 
SPACINGS 


Let A(t) be the area enclosed by the track by time f. So, 
A(t) * SVt. 
Assuming that the area enclosed is approximately circular, 
A(t)=ar'(t). 


Equating these expressions gives 


r(t)* JSVt/z,t>0. (1) 


As illustrated in Figure 5, we can use these expressions to find the track spacing 


(S°), such that r(r) is just tangent to U(t+7)—R at some time (1 ). 


Referring to Figure 5, at time t= 2 
U(t+r)-R=r(t)=JS'Ve/z. 


So, 


=(Ut+(Ur-R)) 


=Ut? +2U (Ur—-R)t+(Ur—R) 


11 


range (nm) 





0 i time, ¢ 


Figure 5. Calculation of S” 


Collecting terms, 


Ur +(2U (Ur—R)-SV/x)t+(Ur-R) =0. 
b c 





This is a quadratic equation in ¢. There can be 0, | or 2 solutions. We are looking 
for the S” that produces one solution. From the quadratic formula, 4ac = b” yields one 


solution. 


So, 


4U?(Ur-R) =(2U(Ur—R)-(SV)/x) 
4U (Ur—R)SV 


=4U?(Ur—R) 
a 


+(SV/z) . 





One solution is $ =0. Assuming that S #0, and then dividing both sides by S, gives 


= = 2 
4nU (Ut ee le 
a a 





Solving for S gives 


12 


se 4aU (Ut—R) 


7 (2) 


S* is the largest track spacing which keeps all search effort on the furthest-on- 


At time t=-—t (initial detection), the target’s location is circular normal with 
means (0, 0) and variances (oer) The target then begins a Random Tour with 
parameters 4 and U (target’s speed). 

By the Central Limit Theorem, the target position at time ft = —z (for sufficiently 
large f) is still circular normal with mean (0, 0), but with variances increased from o (at 
U? (t + T) 

A 


time t=-r) to o° + (at time ¢ >—z ). And since the radial position of a circular 


bivariate normal random variable with variances s” is Rayleigh distributed with mode s, 


we know that the radial position at time f, R(t), is approximately Rayleigh with mode 


oo U Atte) . 
A 


13 


U* (t +T ) distance from datum, r 
A 


o + 


Figure 6. Density Function of Radial Position of Target at Time +, Assuming a Normal 
Approximation to the Random Tour 


So, a spiral search which maintains a radial distance at search time ¢ of 


, U*(r+t) rear ee as 
o° +—————. , t=0 should be a good search. The justification is that this radial 


A 


distance is the mode of the approximate radial distance density function. 


The Archimedes spiral search which provides the tightest support below for this 


function is 


as illustrated in Figure 7. 
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distance 
from origin 










ert (t) = 


oUt, 


time ft 


Figure 7. Calculation of S$ 


Ray 


_ {Ut SVt ! ; 
Setting 7 equal to ,/—— (from Equation (1)) and solving for S yields 
1 


SVt_U’t 
—— =— _,S0 
1 A 
U*x 
a7 - 
Sy Should be a good track spacing for a random tour target that is approximately 


Ray 


normally distributed. However, a small 2 can make S,, unreasonably large because the 


normal approximation may not apply in that case. For instance, when 


10°xz 


s —— —_—— # 15,700 miles, which is 
200 x 0.0001 


U =10kts,V = 200kts, 2 =0.0001/hour , Sey i 


too large. 
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C. S=2R 


When the target is slow or the searcher is fast (i.e., we is small), then an 


exhaustive search (track spacing of 2R) should be effective. 2R is the largest track 


spacing which guarantees no gaps in search coverage. 
D. PROPOSED EXPRESSION 


We have introduced three analytical expressions for the recommended track 


spacing: S S*, and 2R. Each has its weaknesses and strengths. In particular, 


Ray? 


° Spay Should be good when enough time has elapsed to allow the Central 
Limit Theorem to guarantee that the target distribution is roughly bivariate 


normal. But as we have seen, when 4 is small, ae can become 


arbitrarily large. Consequently, it should be bounded above by another 
expression, and we will use S” for this purpose. 


e S” is intended to keep all search effort within the furthest-on-disk. As 


such, S° should be (approximately) the largest track spacing considered. 
Otherwise search effort would be placed in regions where the target can 
not be located. 


e 2R is the largest track spacing which guarantees no gaps in search 
coverage. Thus 2R should be (approximately) the smallest track spacing 
considered for a stationary or slow target. It might, however, be too small 


when R is small or the speed ratio we is large. 


Considering the strengths and weaknesses of S S*, and 2R, the following 


Ray ? 


expression is proposed for the recommended track spacing: 
§ =max{@R,min{BS,,,,5°+7}\. (4) 


Simulation will be used to determine best-fit values of a, (, and vy. The need 
for y was discovered when conducting the simulation experiments. It was noted that as 


the initial observation error standard deviation, o , increased, the best-fit y also increased. 
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V. DESIGN OF EXPERIMENT 


A. NOLH (NEARLY ORTHOGONAL LATIN HYPERCUBE) DESIGN 


Simulation experiments were conducted to find a best-track spacing for the 
Archimedes spiral search pattern. For the design of the experiment, a full factorial design 
for the decision variable S was crossed with a NOLH [5] design for the other five factors. 
We refer to the factors other than S as “Blocking Factors”. They affect the outcomes, and 
we wish to understand their impact, but we have no ability to change them in a real-world 
scenario to achieve a better outcome. To construct the NOLH, a Microsoft Excel 
spreadsheet created by Professor Susan Sanchez (Sanchez, 2005) was used. The number 
of full factorial design points for S is 181 ranging from 0 to 45nm in .25nm increments, 
and the number of NOLH design points is 65. Appendix A is an extract of the NOLH 
design which depicts the 65 design points. The 65 design points are close to orthogonal. 
(All pairwise correlation values for design columns are lower than 0.05.) The two 
designs are crossed to create the following design table. The number of crossed design 


points is 65x 181 = 11,765. 




































































DP U T V R A DP S 
14.656 0.57 153.91 0.992 0.25 1 0 
19.109 | 1.578 | 116.41 | 1.133 | 0.688 

x 2 0.25 
3 | 17.922 | 1.039 | 242.97 | 0.828 | 0.594 
és 3 ans 1 ae T0838 0 = 0 563 a i 
: - - : : Full Factorial 


NOLH Design for Five Blocking Factors Design for Decision Variable S 





DP U T V R A S=0 | S=0.25 es S=45 
1 14.656 | 0.57 | 153.91 | 0.992 | 0.25 
— 2 19.109 | 1.578 | 116.41 | 1.133 | 0.688 
17.922 | 1.039 | 242.97 | 0.828 | 0.594 

















65 3.375 | 1.016 | 109.38] 0.711 | 0.563 
Crossed Design 



































Figure 8. | Crossed Design for Six Factors 
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B. THE MAXIMUM PROBABILITY OF DETECTION (Fd...) FOR EACH 
DESIGN POINT 


Appendix B is the JAVA code for the search simulation model used for data 
collection. One thousand replications are used for each of the 11,765 crossed design 
points. By running the simulation model, the probability of detection (Pd ) for each 


design point was estimated. Table 1 and Figure 9 show some of the Pd results. 








Blocking Pd 
Factor U T V R A 
DP S=0 S=0.25 or S=45 
1 14.656 0.57 153.91 0.992 0.25 0 0.054 i 0.08 
2 19.109 1.578 11641 1.133 0.688 0 0.061 8 0.078 
3 17.922 1.039 242.97 0.828 0.594 0 0.089 ais 0.038 
64 7.531 1.813 144.53 1.156 0.094 0 0.041 ‘ae 0.069 
65 3.375 1.016 109.38 0.711 0.563 0 0.397 8 0.11 


Table 1. Probability of Detection ( Pd ) for Each Blocking Factor Design Point 


Design Point 1 








Pd 


best 











Pd 




















0 10~ best 20 30 40 50 


Figure 9. Plot of Estimated Simulation Probabilities of Detection for Blocking Factor 
Design Point 1 
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Figure 9 shows the simulation Pd for design point 1 and the 181 track spacings 


(0-45nm). We can now determine the largest Pd ( Pd,,.) and the track spacing producing 


best 








Pd yes: (Ses; ) Example results are given in Figure 9 and Table 2. 
Blocking 
Factor U T Vv R A ee PG ox 
DP 
1 14.656 0.57 153.91 0.992 0.25 10.75 0.758 
2 19.109 1.578 116.41 1.133 0.688 22.5 0.116 
3 17.922 1.039 242.97 0.828 0.594 17.75 0.39 
64 7.531 1.813 144.53 1.156 0.094 9.5 0.835 
65 3.375 1.016 109.38 0.711 0.563 1.75 0.943 


Table 2. Table for S,,, and Pd,,., for Each Blocking Factor Design Point 


best best 


C. FINDING THE BEST a, £, AND y 


Using the simulation data of Tables | and 2, we can determine the values of a, 
f, and y in Equation (4) which produce a track spacing and resulting Pd as close as 


possible to the maximum obtainable Pd. We use the following procedure: 


For each of the 51° possible values of (a, B, x); where a, and y vary over 


{0, ple aes es, 5} 
For each of the 65 blocking factor design points (dp) 


Determine $ (a, 8, y, dp) from Equation (4) 


Determine Pd (s ‘ dp) from Table 1 


19 


Determine Pd,,,, (dp) from Table 2 


Determine MOE, (a, 2, y, dp) = P4 hj 


best 


Continue for each of the design points. Then compute 


MOE, (a, Bs: /)=(%]EMOE (a Boys dp) 


Continue for next (a, B, y) to find by enumeration 


(a’, f°, y’) = argmax MOE, (a, f, 7) 


a@.B.7 
For each (a, 8,7) triplet, the analysis generates 65 values of 
MOE, (a@, f, y, dp) — one for each of the blocking factor design points. Then 
MOE, (a, 8, 7) is the average of these 65 values. (a",B",7") are the values which 


maximize the average of the 65 values of MOE,. 
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VI. 


A. MAIN RESULT 


Using the enumeration procedure of the previous section, 51° (= 132,63 1) values 


of MOE, (a, 8, v7) were computed. The largest MOE, was .928, and the maximizing 


arguments were (a’, 8.7’) =(2, 2.7, 1.8). 


SIMULATION RESULTS 


B. VISUALIZATION OF MAIN RESULT USING CONTOUR PLOTS 


Contour plots were constructed to help visualize how MOE, (a, B, 7) varied 


over the range of its arguments. These plots reproduces below indicated that the best 


values of (a, £, y) were approximately (2, 2.7, 1.8). 


Figure 10. 


Contour Plot for Mean Ratio 











beta 














alpha 








Mean Ratio 
MB <_0.200 
MM <_ 0.282 
MB 0.363 
MMB <-0.445 
MM 0.527 
MB <_ 0.608 
ME <_0.690 
MB -0771 
MMB <_ 0.853 
MB <-0.935 
MB. 1.000 


Contour Plot for @ and # given y=0 
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<= 0.216 
<= 0.298 
"<= 0.380 
1 <= 0.461 
MT << 0.543 
( < 0.624 
™ <=0.706 
<= 0.788 
" <£0.869 
7 <0.951 


MB <_ 0.233 
Ml <.0.314 
Ml <- 0.396 
MS <_ 0.478 
lM <. 0.559 
MM. 0.641 
MB... 0.722 
"0.804 
[i <= 0.886 
Ml <. 0.967 


MB <- 0.249 
MB 0.331 
Ml <.0.412 
lS <- 0.494 
lM <.0.576 
MS <_ 0.657 
Ml <. 0.739 
M™ < 0.820 
MS <- 0.902 
MB <. 0.984 


MM <_ 0.265 
MB <. 0.347 
MB <- 0.429 
MB <-0.510 
MM <.0.592 
MB <_ 0.673 
MM <- 0.755 
ll <_ 0.837 
MB <_0.918 
MB <. 1.000 


Figure 11. 


Figure 12. 








Contour Plot for Mean Ratio 
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MM 0282 M2 <o0208 4 MM<o314 MMM<-0.331 MM <-0.347 
MM 0363 = <0380 MM<o396 MMM<-0412 Ml <-0.429 
MMH 0445 © <-0461 M0478 4 M<-0.494 4M <-0.510 
MM 0527 Poss = M<oss9) «96MM 0576 «6M <= 0.592 
© MM <-o6cos %<062 Moe MM<0657 Mil 0.673 
8 MMos © «0706 M0722 MM<o739 60755 
MM <-0771 <0788 3| 1 <0804 Ml <0820 Ml <-0.837 
MM o83 | <ose9 Moses MM<-0902 MM <~o918 
MM <-0935 “%<0951 Ml<o967 M<-0984 Ml <-1.000 
ME. 1.000 
a a a a Oa | T 
0.0 1.0 2.0 3.0 4.0 5.0 
alpha 
Contour Plot for @ and / given y= 1 
_ Contour Plot for Mean Ratio 
Mean Ratio 
MM 0200 MM <o21¢ 8 MM<-0233 MMM<o0249 MMM <-0.265 
MM -.o282 «| <0298 «96MM 0314 «= <= 0.331 «9 MI <= 0.347 
MB <_ 0.363 0.380 MM <.0412 Ml <.0.429 
MB <. 0.445 .461 =0.478 M<o494 Ml <-0.510 
MM -.0527 “8.0543 MM<osso M0576 489M <- 0.502 
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7 MM 0600 3% <o700 9 MM<0722 MM<-0739 Ml <-0.755 
MMeo71 2 <=0.788 «= 0.804 «= <= 0.820 «= MN <-0.837 
MMM <-0.853 | <=0.869 =0.886 MM<o0902 Ml <-0.918 
MMH 0935 38 <0951 M<o967 MMM<o0984 Ml <-1.000 
ME. 1.000 




















Contour Plot for @ and # given y= 2 


22 


Figure 13. 


Figure 14. 





Contour Plot for Mean Ratio 
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Contour Plot for a and # given y=4 
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VII. CONCLUSION / FUTURE RESEARCH 


A. CONCLUSION 


Since the military threat of hostile submarines is increasing, the need for effective 
ASW operations is also increasing. In this thesis, the recommended track spacing for an 
Archimedes spiral search in a datum search problem was studied. This thesis has sought a 
meta-model to recommend a good track spacing. This analysis combined three analytical 


functions into a single parameterized expression. To find the best-fit parameters 


maximizing the average ratio BP yy, d 8 simulation experiment with a NOLH (Nearly 
best 
Orthogonal Latin Hypercube) design was used. The complete model is 


§ =max{2R, min{2.7S_,S° +1.8}}, 


ray? 


where R is detection range, S° and S,,. are given by Equations (2) and (3). 


Ray 
B. FUTURE RESEARCH 


This thesis blended different techniques involving stochastic modeling and 
simulation, design of experiments, and optimization to create a meta-model for track 


spacing for an Archimedes Spiral. Suggested future research includes: 


e Good search patterns with multiple targets and multiple searchers; 

e Selection of appropriate spiral search pattern in datum search problem; 
e Factor analysis for probability of detection in ASW datum search; 

e Recommended track spacing when the other parameters are unknown; 
e Impact on observation error for probability of detection. 
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APPENDIX A. NOLH (NEARLY ORTHOGONAL LATIN 
HYPERCUBE) DESIGN FOR FIVE BLOCKING FACTORS 


Low level 1 0.5 100 0.5 0 
High level 20 2 250 2, 2 
decimals 3 3 24 3 3 
factor name U T Vv R A 
1 14.656 0.57 153.91 0.992 0.25 
2 19.109 1.578 116.41 1.133 0.688 
3 17.922 1.039 242.97 0.828 0.594 
4 13.172 1.836 207.81 1.18 0.125 
5 18.516 1.203 128.13 0.523 0.188 
6 11.094 1.883 135.16 WA 0.313 
7 15.547 0.781 179.69 0.547 0.5 
8 16.438 1.648 238.28 0.969 0.75 
9 14.063 0.547 102.34 1.719 0.813 
10 19.406 1.531 172.66 1.648 0 
11 10.797 0.523 245.31 1.297 0.406 
12 19.703 1273 203.13 1.883 0.344 
13 11.391 0.828 149.22 1.461 0.781 
14 15.844 32 163.28 1.836 0.531 
15 11.984 0.992 217.19 1.93 0.875 
16 13.766 1.367 182.03 1.578 0.281 
17 17.328 1.156 158.59 0.688 1.344 
18 11.688 1.789 139.84 1.016 1.156 
19 14.953 1.086 184.38 0.945 1.563 
20 12.578 1.602 228.91 0.594 1.094 
21 16.734 0.805 160.94 0.641 1.969 
22 16.141 1.742 100 0.734 1.031 
23 20 1.109 219.53 0.875 1.938 
24 12.281 2 Diy ISD 1.063 1.375 
25 12.875 0.875 114.06 1.391 1.625 
26 15.25 1.438 123.44 2 1.531 
27 17.031 0.641 2125 1.695 1.063 
28 18.219 1.766 193.75 1.742 1.844 
29 18.813 0.945 118.75 1.414 1.281 
30 14.359 1.906 151.56 1.602 1.781 
31 13.469 0.688 205.47 1.344 1.906 
32 17.625 1.484 240.63 1.789 1.438 
33 10.5 N25 175 1.25 1 
34 6.344 1.93 196.09 1.508 eS 
35 1.891 0.922 233.59 1.367 1.313 
36 3.078 1.461 107.03 1.672 1.406 
37 7.828 0.664 142.19 1.32 1.875 
38 2.484 1.297 221.88 1.977 1.813 
39 9.906 0.617 214.84 1273 1.688 
40 5.453 1.719 170.31 1.953 15 


Low level 1 0.5 100 0.5 0 
High level 20 2 2 2 
decimals a 3 2 3 3 
factor name © U t R a 
41 

42 

43 

44 

45 

46 
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APPENDIX B. JAVA CODE FOR DATA COLLECTION 


The Java code within this appendix collects all the simulation results over 11765 


crossed design points. 
/ * 


* File: Arch.java 
* Created on Sep 15, 2007, 12:00 PM 
ay 


[** 
* @author Byungsoo Son 
* Class to collect the probability of detection 
* for the 11765 design points 
*/ 


import java.util.Random; 
public class Arch { 


public static void main(String[] args) { 
//set input values 
Random generator = new Random(); 
//initializing all variables 
double targetSpeed=0; 
double tau=0; 
double tMax =5; 
double searcherSpeed=0; 
double detectionRange=0; 
double trackSpacing=0; 
double lambda=0; 
double timeStep=0; 
double xTarget=0; 
double yTarget=0; 
double sigma=0; 
double xSearcher=0; 
double ySearcher=0; 
double distance=0; 
int Detection=0; 
double courseDirection; 
int numberRelications = Integer.parseInt(args[0]); 
int numberDetection=0; 








if (args.length == 9) { 
sigma = Double.parseDouble(args[1]); 
targetSpeed = Double.parseDouble(args[2]); 
tau = Double.parseDouble(args[3]); 
searcherSpeed = Double.parseDouble(args[4 
detectionRange = Double.parseDouble (args [ 
lambda = Double.parseDouble(args[6]); 
trackSpacing = Double.parseDouble(args[7]); 
timeStep = Double.parseDouble(args[8]); 











1); 
5) )% 

















} 
else { 
System.err.println("Supply nine arguments: " + 
"numberReplications(int), sigma(double) " + 
"targetSpeed(double), tau(double), " + 
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"searcherSpeed(double), " + 
"detectionRange (double), lambda(double), " + 
"trackSpacing(double) " + 
"and timeStep(double)"); 
System.exit(-1); 
} 
// replicate n times given the other parameters 
for(int n=0;n<numberReplications;nt+t) { 
xTarget = generator.nextGaussian() *sigma; 
yTarget= generator.nextGaussian() *sigma; 
xSearcher=0; 
ySearcher=0; 
Detection=0; 
courseDirection = 
generator.nextDouble() *2*Math.PI; 
//target motion model for time -tau to 0 
for(int i=0;i<(tau/timeStep) ;i++) { 
if (generator.nextDouble()<lambda*timeStep) { 
courseDirection = 
generator.nextDouble()*2*Math.PI; 








} 








xTarget = xTarget + Math.cos(courseDirection) 
*xtargetSpeed*timeStep; 
yTarget = yTarget + Math.sin(courseDirection) 


*targetSpeed*timeStep; 





} 
//searcher and target motion model for time 0 to tMax 
//T£ the distance between searcher and target is less 
//than the detection range of sensor,detection occurs 
for(int j3=0; 3<(tMax/timeStep) ; j++) { 
if (generator.nextDouble()< lambda*timeStep) { 
courseDirection = generator.nextDouble() 
*2*Math.PI; 
} 
xTarget = xTarget + Math.cos(courseDirection) 
* targetSpeed*timeStep; 
yTarget + Math.sin(courseDirection) 
* targetSpeed*timeStep; 
Math.sqgrt(trackSpacing * searcherSpeed 
* (timeStep * J) /Math.PTI) 
* Math.cos(2*Math. sqrt (searcherSpeed 
*Math.PI* timeStep*j/trackSpacing) ); 
Math.sgrt(trackSpacing * searcherSpeed 
* (timeStep * j)/Math. PTI) 
* Math.sin(2*Math. sqrt (searcherSpeed 
*Math.PI * timeStep*j/trackSpacing) ); 
distance = Math.sqrt (Math. pow(xTarget—xSearcher, 2) 
+ Math.pow(yTarget-—ySearcher,2)); 
if(distance < detectionRange) { 
Detection = Detection + 1; 
break; 

















ylarget 





xSearcher 


ySearcher 








} 
} 


numberDetection = numberDetection + Detection; 





} 


System. out.println( (double) numberDetection/numberReplications) ; 
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